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ABSTRACT 

We present the discovery of three new Milky Way satellites from our search for compact stellar 
overdensities in the photometric catalog of the Panoramic Survey Telescope and Rapid Response 
System 1 (Pan-STARRS 1, or PSI) Stt survey. The first satellite, Laevens 3, is located at a heliocentric 
distance of d = 67 ± 3kpc. With a total magnitude of My = —4.4 ± 0.3 and a half-light radius 
Th = 7 ± 2pc, its properties resemble those of outer halo globular clusters. The second system, 

Draco H/Laevens 4 (Dra H), is a closer and fainter satellite (d ^ 20kpc, My = —2.9 ± 0.8), whose 
uncertain size (r/j = 191® pc) renders its classification difficult without kinematic information; it 
could either be a faint and extended globular cluster or a faint and compact dwarf galaxy. The 
third satellite, Sagittarius H/Laevens 5 (Sgr H), has an ambiguous nature as it is either the most 
compact dwarf galaxy or the most extended globular cluster in its luminosity range (r^ = 37lg pc and 
My = —5.2 ± 0.4). At a heliocentric distance of 67 ± 5kpc, this satellite lies intriguingly close to the 
expected location of the trailing arm of the Sagittarius stellar stream behind the Sagittarius dwarf 
spheroidal galaxy (Sgr dSph). If confirmed through spectroscopic follow up, this connection would 
locate this part of the trailing arm of the Sagittarius stellar stream that has so far gone undetected. It 
would further suggest that Sgr H was brought into the Milky Way halo as a satellite of the Sgr dSph. 

Subject headings: Local Group — Milky Way, satellites, dwarf galaxies, globular clusters, streams: 
individual: Sagittarius H, Draco H, Laevens 3 


1. INTRODUCTION 

Two decades ago the prevalent view of the Milky 
Way (MW) as an isolated system was radically changed 
by the discovery of a tidally disrupting dwarf galaxy, 
embed ded in a stream in the constellation of Sagit¬ 
tarius (|Ibata et al.l II994I1 . highlighting the underrated 
importance of Milky Way-satellite interactions. With 
ACDM models predicting a whole new population of 
faint satellite dwarf galaxies (DCs) orbiting the MW 
(e.g. iBullock et al.l[2000l 1200 Ih . the new challenge was 
to find these, until then, elusive objects. At the turn 
of the century, the advent of large CCD surveys such as 
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the Sloan Digital Sky Survey (SDSS; lYork et Ml 12000(1 
uncovered some 16 dwarf galaxies among the faintest 
ever found fe.g. IWillman et al.ll2005t iZucker et al.ll2006l: 


iBelokurov et alJ 120071 : iWalsh et al.l 1200711 . Though rev¬ 
olutionizing our view of the satellite galaxies, just a 
handful of n ew globular clusters (GCs) were found, faint 
and nearby (Kodosov et al.ll2007t fBelokurov et al.ll201fll : 


iMuhoz et al. 120121 : iBalbinot et al. I2013f) . In addition, 
the SPSS enabled the discovery of several tidal strearn s 
(e.g. IBelokurov et al.l[2006l : iGrillmair fc Dionatosll2006ll . 
further illustrating the central role of satellite and cluster 
disruption in building up the MW’s halo. 

With the second generation of surveys emerg- 
ing such as the Dark Energy S urvey (DES, 
iThe Dark Energy Survey Collaborationl l2005f) . PSI 
(Chambers et al. in preparation), later data releases of 
the SDSS, the Survey of the MAgellanic Survey History 
(SMASH, Nidever et al. in preparation), and VST Atlas 
have seen the number of known MW likely dwarf galaxies 
expand further from ^25 to ^35. These once elusive 
systems appear to be more commo n as deeper, but also 


wider coverage, data are gathered (iLaevens et ahlPOTa 
Bechtol et al.l [201^ iKopo sov et al ] 20151: Martin et al l 


20151 : IKim et al.l I2015al : IKim fc Jerienl 120151^ In par- 

allel, these systematic surveys als o revealed a smaller 
number of faint n ew GCs (e.g. ILaevens et al.l 12014 
IKim &: Jerienl 1201511 . The increase in the number of 
MW satellites led to the blurring of the traditional 
distinction between small, baryon-dominated GCs and 
more extended, dark-matter dominated DCs. Taking 
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the photometric properties of these new satellites at 
face value shows that they straggle the DG and GC 
boundary in the size-lu minosity plane, in th e so-called 
“valley of ambiguity” (jGilmore et al.l l20(T^ . Though 
follow-up observations have implied velocity dispersions 
higher than that expected from the stellar content 
for most of the ii e w satellites (iMartin et al.l 120071: 
Simon fc Gehal l2007t iWillman et al.l 120111 : 1s imon et alJ 
201 ll: iKirbv et al.l 1201311 . those measurements suffer 
from small number statistics and the unknown effect of 
binary stars on the kinema tics of these small systems 
(jMcGonnachie fc G6t^l2f)Tol l. 

The recent discoveries of such faint candidate DGs out 
to ~ 70kpc within DES confirm that they are in fact 
common and that they could indeed correspond to the 
large population of faint dark-matter dom inated systems 
expected to inhabit the M W halo (e.g. iTollerud et alJ 
[2?iM IBullock et ai] [Mflil . These new satellite s, lo¬ 
cated close to the Magellanic Glouds (iBechtol et alJlMT^ 
IKirbv et al.ll2015t iKoposov et alJl2015ll . have emphasized 

the tendency of these faint stellar systems to be brought 
into the MW surroundings in groups. Moreover, appar- 


entlv isolated systems often share a 

Droximitv with stellar 

streams (e.g. 

Belokurov et al.l 12008 

: IDeason et al.l 12014 

ILaevens et al. 

20151: IMartin et al.l 201511. 


Over the last three years, our group has focused on 
the search for compact stellar systems in PSl, which 
has so far revealed two new MW satell i tes: a likely GC, 
Laev ens 1/Crater (jLaevens et al.l[20H: iBelokurov et al.l 
[2?fn. as well as a very faint satellite Triangulum II 
(jLaevens et al.l I20I5I1 , whose nature has not yet been 
confirmed by spectroscopy. In this paper, we present the 
discovery of three new MW satellites discovered from the 
latest PS I photometric catalog: a faint GC, Laevens 3 
(Lae 3); a faint satellite, Draco II/Laevens 4 (Dra II), 
whose uncertain properties make its nature ambiguous; 
and another ambiguous system, Sagittarius II/Laevens 5 
(Sgr iijElB. This paper is structured in the following 
way: in section 2 we describe the PSl survey and briefly 
outline the method which led to the discovery of the three 
satellites. In section 3 we discuss the properties of Lae 3, 
Dra II, and Sgr II, concluding and discussing the impli¬ 
cations of the discoveries in section 4. 

In this pape r , all in agnitudes are dereddened using the 
iSchlegel et al.l (Il998fl maos, adopting th e extinction co- 
efScients of iSchlaflv &: Finkbeined (|2Cillll . A heliocentric 
distance of 8 kpc to the Galactic center is assumed. 

2. THE 37r PSl SURVEY AND DISCOVERY 

The PSl survey (K. Chambers et al. in preparation) 
observed the whole sky visible from Hawaii {6 > —30°), 
providing an unparalleled panoptic view of the MW and 
the Local Group. Throughout the 3.5 years of the Stt 
survey, the 1.8 m PSl telescope, equipped with its 1.4- 
gigapixel camera capable of observing a 3.3-degree field of 
view, collected up to four exposures p er year in five differ - 
ent optical filters: fom tpApi zpi r/pi: iTonrv et al.l[2M^ . 

We assign double names for these last two systems, pending 
spectroscopic confirmation as to the nature of these stellar systems. 
For convenience and clarity, throughout the remainder of the paper, 
we refer to these satellites by their constellation name 

We assign Roman numeral II to this system and refer to 
the Sa gittarius dwarf spheroidal galaxy discovered bv llbata et al.l 
l|1994l~) as Sgr dSph. 


Once the individual frames have been taken at Haleakala 
and downloaded from the summit, the photomet ry is dea¬ 
erated through the Ima^e Proce ssing; Pipeline (jMagnieri 
[2?M[2?i(T7l:fMa,gnier el RLl[2nnl . 

The internal 37r stacked catalogs were released in three 
processing versions (PV), with each consecutive version 
corresponding to a higher number of individual exposures 
and improved photometry. The three stellar systems de¬ 
scribed in this paper were found using the intermediate 
PV2 catalog and supplemented with the upcoming PV3 
photometry for their analysis, when beneficial. Although 
there are many small differences between the two pro¬ 
cessing versions, their most interesting features for our 
study are that the PV2 psf photometry is performed on 
the stacked images, whereas for PV3, the stacks are only 
used to locate sources before performing the photometry 
on each individual sub-exposure, with its appropriate psf. 
As a consequence, the PV3 photometry is more accurate, 
but the PV2 star/galaxy separation is more reliable. The 
depths of the bands of PV2, enabling the discoveries, are 
comparable to the SDSS for the gn band (23.0) and 
reach ~ 0 5/^ 1.0 magnitude deeper for rpi (22.8) and 
ipi 122.5: [Metcalfe et'ffll2013[ l. 

With large GGD surveys, automated search algorithms 
were developed to perform fast and efficient searches 
of these massive data sets for the small stellar over¬ 
densities that betray the presence of faint MW satel¬ 
lites. Thes e techniques, originally implemented o n the 
SDSS data (IKoposov et al.l[2n?iM iWalsh et al.ll2009l l have 
proven very successful. Inspired by this, we have devel¬ 
oped our own similar convolution technique (Laevens et 
al. in preparation), adapted to the intricacies of the 
PSl survey. The technique consists in isolating typi¬ 
cal old, metal-poor DG or GC stars using de-reddened 
color-magnitude information [(r — z)o,fo]- For a chosen 
distance, masks in color-magnitude space are determined 
based on a set of old and metal-poor isochrones. The dis¬ 
tribution of sources thereby extracted from the PSl stel¬ 
lar catalog is convolved with two different window func - 
tions or Gaussian spatial filters (jKoposov et al.l l^008jl . 
The first Gaussian is tailored to the typical dispersion 
size of DGs or GCs (2', 4', or 8'), whereas the second 
one accounts for the slowly-varying contamination on far 
larger scales (28' and 56'). Subtracting the map pro¬ 
duced from convolving the data with the larger Gaussian 
from that obtained with the smaller Gaussian results in 
maps of the PSl sky tracking over- and under-densities 
once we further account for the small spatial inhomo¬ 
geneities present in the survey. After cycling through dif¬ 
ferent distances and the aforementioned sizes, we convert 
and combine all the density maps into statistical signifi¬ 
cance maps, allowing for a closer inspection of highly sig¬ 
nificant detections that do not cross-match with known 
astronomical objects (Local Group satellites, background 
galaxies and their GC systems, or artifacts produced by 
bright foreground stars). We further weed out spurious 
detections, by checking that these over densities do not 
corr espond to significant background galaxy overdensi¬ 
ties (IKoposov et al.ll2(i0^ . Applied to PVl, this method 
already led to the discovery of th e most distant MW 
globular cluster Laev ens 1/Crater (jLaevens et al.ll20l^ 
IBelokurov et al.ir2014ll . a s well as one of the f aintest MW 
satellites, Triangulum II (iLaevens et aLll^lSH . whose na¬ 
ture is not yet known. Sgr II, Dra II, and Lae 3 were 
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detected as 11.9, 7.4, and 6.5 a detections, comfortably 
above our 5cr thresholcQ- All three new satellites lie out¬ 
side the SDSS footprint, which explains why they were 
not discovered before. Sgr II and Lae 3 are located at 
fairly low Galactic latitudJ*^ (6 ~ —20°) and Dra II is 
quite far north {6 ~ -1-65°). 

3. PROPERTIES OF THE THREE STELLAR SYSTEMS 

3.1. Color-Magnitude Diagrams and Distances 
3.1.1. Sagittarius II 

The color-magnitude diagram (CMD) of stars within 
one half-light radius of Sgr II (2.0'; see below for the 
structural parameters) is displayed in the top row of pan¬ 
els in Figure [U next to the CMD of a field region of the 
same size. Since the features of Sgr II are so obvious 
we rely here only on the more accurate PV3 photom¬ 
etry at the cost of a poorer star/galaxy separation at 
the faint end. Given the location near the MW bulge 
[(£,6) = (18.9°,-22.9°)] the field CMD is very pop¬ 
ulated, but the Sgr II CMD features are nevertheless 
clearly defined with a red giant branch (RGB) visible 
between [(rpi —ipi)o,ipi,o] — [0.30,16.5] and [0.15,21.5], 
before its main sequence turnoff at Jpi.o < 22.0. The 
most obvious feature, however, is the horizontal branch 
(HB) of the system, clearly visible for 19.5 < ipi_o < 20.0 
and (rpi — ipi)o < 0.0. When selected with the box over¬ 
laid in orange on the CMD of the second Sgr II CMD 
panel of FigurelU these stars correspond to a well-defined 
spatial overdensity (right-most panel). Isolating the HB 
stars (blue box) also highlights how clustered they are on 
the sky. We further highlight a single star that is bluer 
and brighter than the turn off and could potentially cor¬ 
respond to a blue straggler (green box and green trian¬ 
gle). 

The presence of the reasonably well-populated HB, 
with 13 stars within 3 half-light radii, allows for a ro¬ 
bust esti mation of the distan ce to the satellite. Equa¬ 
tion 7 in (jPeason et al.ll20lTll describes the relation be¬ 
tween the absolute magnitude of these HB stars and 
their SDSS g — r colors. Converting the PSl magni¬ 
tudes to the SDSS bands, for which the relation holds, 
reveals a median g = 19.60 ± 0.03 and Mg = 0.47 ± 0.04 
when we perform a Monte Carlo resampling of the 
stars’ uncertainties. These lead to a distance-modulus 
of 19.1 3 ± 0.15, wher e an u ncertainty of 0.1 was assumed 
on the lDeason et al.l (|2014| 1 relation. This translates into 
a Heliocentric distance of 67 ± 5kpc or a Galactocentric 
distance of 60 ± 5 kpc. Fixing the satellite at this dis¬ 
tance modulus, we experiment with isochrones. Overlaid 
on the Sgr 11 CMD of Figure [U we also show the old and 
metal-poor isoch rone from the Pars ec library (12Gyr, 
[Fe/H] = —2.2; iBressan et al.l[2?)T^ that provides the 
best qualitative fit to the CMD features at this distance. 

The properties of Sgr 11 are summarized in Table 1. 

3.1.2. Draco II 

Draco H is much closer and less luminous than Sgr H, 
as can be seen in the CMD of stars within 2rh of the 

For context, applying this technique leads to the recovery of 
the faint satellites Segue 1 and Bootes I, originally discovered in 
the SDSS, with a significance of 7.3 and ll.hcr, respectively. 

In fact, Lae 3 is clearly visible on the DSS plates and could 
have been discovered before the PSl era. 


satellite’s center in the second row of panels in Figure [TJ 
Here, since we need both depth and a good star/galaxy 
separation to clean the main sequence of Dra H, we use 
the PV3 photometry combined with the superior PV2 
star/galaxy flags. This has the consequence of remov¬ 
ing some faint PV3 stars misidentified as galaxies but 
more optimally cleans the main sequence of the satel¬ 
lite. A field CMD is shown in the right-most CMD panel 
and helps identify the Dra H features: a populated main 
sequence between [(rpi — *pi)o, *Pi,o] — [0.0,20.2] and 
[0.2,22.0]. At brighter magnitudes, Dra H shows no HB 
and no prominent RGB. However, a group of stars at 
[0.2,17.0] is compatible with being the system’s sparsely 
sampled RGB. As for Sgr H, isolating the stars in these 
CMD features (orange box in the central CMD panel) 
highlights the stellar overdensity in the spatial distribu¬ 
tion shown in the right-most panel. As for Sgr H, we 
identify a potential blue straggler in g reen. 

Due to the absence of any HB staj^T we cannot reliably 
break the distance-age-metallicity degeneracy with the 
PSl data alone. Consequently, we explored isochrones 
of different ages and metallicities, located at varying dis¬ 
tance. The best fit is provided by the Parsec isochrone 
shown in Figure [T] it has an age of 12Gyr and [Fe/H] = 
—2.2 and is located at a distance modulus of 16.9 ± 0.3 
but we caution the reader on the reliability of this par¬ 
ticular isochrone that needs to be confirmed from deeper 
data. 

3.1.3. Laevens 3 

As can be seen in Figure [H Lae 3 is a compact stel¬ 
lar svsterrP^. As such, the automated PSl pipeline fails 
to extract the photometric information of the central 
region that suffers from crowding. We therefore per¬ 
form custom photometry of t his sky cell using dao phot, 
using the same method as (iLaevens et al.ll2014fl . The 
resulting CMD for stars within 3rh of the stellar sys¬ 
tem’s centroid is shown on the bottom row of panels 
in Figure [T] This CMD is still likely to suffer from 
crowding, yet it reveals features that are clearly not ex¬ 
pected in the field population: the Lae 3 RGB between 
[(^■pi — *pi)o)*Pi,o] — [0.75,18.0] and [0.15,21.0], followed 
by the system’s MSTO at fainter magnitudes. The stars 
between [(rpi —ipi)o, *pi,o] — [0.0,16.2] and [0.2,18.2] are 
foreground contaminants and are situated far from the 
satellite center, just under the Mh limit. As for the two 
other satellites, selecting these stars only (orange box in 
the middle CMD) highlight a clear stellar overdensity. 

An investigation into the presence of RR Lyra stars in 
the PSl temporal data (Hernitschek et al. 2015, in prep.; 
Sesar et al. 2015, in prep.) reveals one obvious can¬ 
didate, 0.6 arcmins away from the center of the cluster 
(highlighted by the blue box in the middle CMD and rep¬ 
resented by a blue star in the right-most panel). Briefly, 
RR Lyrae stars are identified in PSl data by providing 
average PSl colors and various va riability statistics to a 
trained Random Forest classifier ([Richards et al.l [201111 . 
The resulting RR Lyrae sample is 80% complete (up to 

This is not per se surprising as, for instanc e, the similarly 
faint system Willman 1 only contains two HB stars llWillman et al.l 

HoH. 

Note that the PSl postage stamp images show no clear stellar 
over density for Sgr II and Dra II, hence why we do not include 
them. 
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Figure 1. Left: from top to bottom, CMD of stars within 1 half-light radius (Sgr II), 2 half-light radii (Dra II) or 3 half-light radii (Lae 3) 
with the favored isochrone: 12 Gyr and [Fe/H] —2.2 for Sgr II and Dra II, and 9 Gyr and [Fe/H] = —1.9 (Lae 3). In the case of Lae 3, we 
also show the same old and metal—poor isochrone (green) as for Sgr II and Dra II, highlighting that with current photometry, it is difficult 
to distinguish between the two. Middle-left: Same as left-most panel with the CMD selection box used to isolate the RGB, HB, MSTO, 
and/or MS stars of the satellites (orange), and an HB selection box (blue) for Sgr 11. The RR Lyra star for Lae 3 is highlighted in blue. 
Candidate blue straggler stars are identified for Sgr II and Dra II (green). Middle-right: CMD of field regions for stars 15 arcmin West 
of the satellites of similar sizes to those used for the left-most panels. Right: Spatial distribution of all stars around the three satellites 
(light black dots) and of stars selected with the orange CMD selection boxes in the middle-left panel (big black dots). For Sgr II, HB stars 
corresponding to the blue HB selection box in the middle panel are represented by blue stars symbols. Finally, candidate blue straggler 
stars corresponding to the green box are displayed as green triangles. 
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Figure 2. : gpirpiipi image of Laevens 3, built from the stacked 
PV3 PSl images. The image is 2.5' x 2.5' and North is to the top 
and East is to the left. 


80 kpc) and 90% pure. The distances of PSl RR Lyrae 
stars are uncertain at the 5% level. The RR Lyra star 
in Lae 3 has also been observed more t han 100 times by 
the Palomar Transient Factory (PTF; iLaw et ^ 120091 : 
iRau et all 1200911 . The distance and period measured 
from PSl data agree within 2 kpc and 5% with those 
measured from the PTF data. The RR Lyra star is at 
(rn — M)o = 19.14±0.10, or a distance of 67±3kpc. The 
~14 hour period of the star suggests a star with a metal- 
licity range of —1.9 < [Fe/H] < —1.6. Fixing this dis¬ 
tance, we once again experiment with various isochrones 
and conclude that the CMD feature are best tracks by an 
isochrone with a comparatively young age of 9 Gyr and 
[Fe/H] ^ —1.9, compatible with the properties of the RR 
Lyra star. The isochrone fit, fixed at that distance, tracks 
the main features of the satellite such as the RGB and 
the MSTO. Though two blue stars are also present at the 
same magnitude as the tentative red HB, the bluest of 
the two is a field variable star, incompatible for being a 
member of Lae 3. We further comp are the CMD features 
with the GC fiducial published bv iBernard et al.l ()2ni4r i 
in the PSl photometric system. The fiducials of globu¬ 
lar clusters NGC 1904, NGC 5897 and NGC 7089, with 
— 1.9 < [Fe/H] < —1.6 provide a good fit to the CMD 
features and confirm our impression from the isochrones. 

The derived properties are summarized in Table 1 for 
all three satellites. 

3.2. Structural Parameters and Absolute Magnitudes 

The structural parameters of Sgr 11, Dra 11, and Lae 3 
are derived using a version of the iMartin et al.l (l2008f l 
likelihood technique updated to a full Markov Chain 
Monte Carlo framework (Martin et al. 2015, in prep.). 
Using the star’s location in the vicinity of satellite, the al¬ 
gorithm calculates the posterior probability distribution 
function (PDF) of a family of exponential radial density 
profiles, allowing for flattening and a constant contam¬ 
ination from field stars. The parameters of the models 
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Figure 3. Left: Probability distribution functions for the elliptic- 
ity (e), the position angle (9), and the angular and physical half- 
light radii, of Sgr II (top), Dra II (middle), and Lae 3 (bottom). 
Right: Comparison between the favored radial distribution profile 
(full line) and the data, binned according to the preferred struc¬ 
tural parameters (dots), selected as the mode of the PDFs (grey 
line in the left-hand panels). The error bars assume Poissonian 
uncertainties and the dashed line represents the field density. 


are the centroid of the system, the ellipticity, m, the 
position angle, 9 (defined as the angle of the major axis 
East from North), the half-light radius r^O, and finally 
the number of stars, N* within the chosen CMD selec¬ 
tion box. We further determine the physical half-light 
radius from the angular one by randomly drawing dis- 


The ellipticity, here, is defined as e = 1 — b/a with a and b 
the major and minor axis scale lengths, respectively. 

In fact, the algorithm constrains the half-density radius, but 
this is similar to the more common half-light radius if there is no 
mass segregation in the system. 
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Table 1 

Properties of Laevens 3, Draco II and Sagittarius II 



Laevens 3 

Draco II/Laevens 4 

Sagittarius II/Laevens 5 

a (ICRS) 

21:06:54.3 

15:52:47.6 

19:52:40.5 

5 (ICRS) 

+14:58:48 

+64:33:55 

-22:04:05 

en 

63.6 

98.3 

18.9 

bn 

-21.2 

+42.9 

-22.9 

Distance Modulus 

19.14 + 0.10 

~ 16.9 + 0.3 

19.13 + 0.15 

Heliocentric Distance (kpc) 

67 + 3 

20 + 3 

67 + 5 

Galactocentric Distance (kpc) 

64 + 3 

22 + 3 

60 + 5 

My 

-4.4+ 0.3 

-2.9+ 0.8 

-5.2 + 0.4 

Ly 

j^q3.7±0.1 

j^q3.1±0.4 

04.0+0.1 

[Fe/H] 

- 1.9 

~ -2.2 

~ -2.2 

Age (Gyr) 

ro 9 

~ 12 

~ 12 

E{B - y)“ 

0.073 

0.016 

0.097 

Ellipticity 

0.21 + 0.21 

24 

u.zo_Q 23 

Position angle (from N to E °) 


70 + 28 

79+28 

'^-20 

(arcmin) 

0-40t°;°I 

9 7+1-0 
^•'-0.8 

^•'-'-0.3 

rh (pc) 

7 + 2 


38):® 


“ from ^chlegel_e^^l^ and ^chl^^^^inkbeiner (201^ 


tances from the distance modulus values. 

The PDF for the ellipticity, position angle, as well as 
the angular and physical half-light radii may be seen in 
Figure [3] for, from top to bottom, Sgr II, Dra II, and 
Lae 3. The three right-most panels of the figure com¬ 
pare the radial profile of a given satellite, binned follow¬ 
ing the favored centroid, ellipticity, and position angle to 
the favored exponential radial density profile; they dis¬ 
play the good quality of the fit in all cases. All three 
systems are rather compact, with angular half-light radii 
of 2.0^0 3 , 2.7lj!) g, and 0.401 q'°]^ arcmin for Sgr II, Dra II, 
and Lae 3, respectively. However, the different distances 
to these systems lead to different physical sizes: 38l®, 
19^6, and 7±2 pc. In all three cases, the systems appear 
mildly elliptical but the PDFs show that this parameter 
is poorly constrained from the current data. It should be 
noted that, in the case of Lae 3, the crowding at the cen¬ 
ter of the stellar system could lead to an underestimation 
of the compactness and luminosity of the system. How¬ 
ever, the Lae 3 radial profile shows no sign of a central 
dip. 

The absolute magnitude of the three stellar sys¬ 
tems was determined using the same procedure as fo r 
Laevens I and Triangulum H (jLaevens et al.l[^ l4l2015h. 
as wa s also described for the first time in iMartin et al.l 
(|2008ll . Using the favored isochrones and their associ¬ 
ated luminosity functions for the three satellites, shifted 
to their favored distances, we build CMD pdfs after fold¬ 
ing in the photometric uncertainties. Such CMDs are 
populated until the number of stars in the CMD selec¬ 
tion box equals the favored number of stars N* as deter¬ 
mined by the structural oarameterJ*^. The flux of these 
stars is summed up, yielding an absolute magnitude. In 
practice, this operation is repeated a hundred times with 

For this part of the analysis, we make sure to use a selection 
box that remains ~ 1 magnitude brighter than the photometric 
depth so the data is close to being complete. In the case of Lae 3, 
the sparsely populated CMD prevents us from doing so as we need 
to use the full extent of the CMD to reach convergence in the struc¬ 
tural parameter analysis. This is likely to slightly underestimate 
the luminosity of the cluster. 


different drawings of the Markov chains, thus taking into 
account three sources of uncertainty: the distance mod¬ 
ulus uncertainty, the uncertainty on the number N* of 
stars in the CMD selection box, and shot-noise uncer¬ 
tainty, originating from the random nature of populat¬ 
ing the CMD. This procedure yields total magnitudes in 
the PSl rpi band, which we then convert to the more 
commonly used U-band magnitudes through a constant 
color offset {V — r = 0.2) determined from the analy¬ 
sis of more populated, known, old and metal-poor MW 
satellites. This yield My = —5.2 ± 0.4, —2.9 ± 0.8, and 
—4.4 ± 0.3 for Sgr H, Draco H, and Laevens 3, respec¬ 
tively. All three systems are rather faint, as expected 
from their sparsely populated CMDs. 

4. DISCUSSION 

Figure |4] displays the properties of the three new 
discoveries in the context of the other MW satellites 
(GCs or DCs). The top panel shows the size-luminosity 
plane while the bottom panel focuses on the distance- 
luminosity plane. These parameters can already be a 
first indicator as to the nature of these objects, which 
we proceed to discuss here as well as the possible stream 
associations these objects may have. 

4.1. Sagittarius II 

Sgr H occupies an interesting place in the vs. My 
plane as it lies in the very middle of th e “val ley of 
ambigu itj0’ highlighted by iGilmore et al.l ()2007[1 . Al¬ 
though other MW satellites are known with similar ab¬ 
solute magnitudes, Sgr H is smal l er tha n Goma Berenices 
(rh = 74 ± 4 pc: iMunoz et al.l I20I0L Pisces H (rh = 
58 ± I0pc:lSand et al.ll20I2tl . Hvdra H (ru = 68 ± 11 pc; 
IMartin et al.ll201^ or the even larger Leo IV and Le o V 
(rh = 205 ± 36 and 133 ± 31 Dc: lde Jong et al.ll20ir)fl . or 
Ursa Major I {rh = 318 ± 45 pc; IMartin et al.ll200^ On 
the other hand, Sgr H re mains larger than the largest 
GG, Pal 14 (r/i ^ 25 pc; iHilkeil [2006f) . or the recently 

The region in vs. My space that straggles the ‘classical’ 
boundaries between DGs and GCs. 


























Three new Milky Way satellites in PSl 


7 


o 




LJ 


00 

d 


(O 

d 


CNJ 

d 


o 


rh (pc) 

o 




LJ 


00 

d 


(£> 

d 


d 


CM 

d 


o 


Heliocentric Distance (kpc) 

Fi gure 4. Top: Distribution of MW satellites in the size-luminosity plane, color- coded by their elliptici ty. Squares represent GCs from 
the IHarrli II2010I) catalog, supplemen ted by the more recent discoveries of Segue 3 JBelokurov et aril2010l) . Munoz 1 d Munoz et al .11^0121) . 
and Balbinot 1 dBalbinot et al1l2Q131 . Milk y Way confirm e d dwa rf ga laxies are shown as ci rcled dots, with their properties taken from 
IMcConnachid l|2012l '). The co-discoveries by IBechtol et alJ l|2015h and IKodosov et al.l l|2015l^ are shown with triangles and filled circles 
respectively, with t he co-discoveries l inked to each other by a black solid line reflecting the two groups’ different measurements. The 
I Kim et al.l ll2Q15alb IKim et al.l d2Q15bl' ).and IKim &; JerienI d 20 151) satellites are shown with diamonds. Hydra II, discovered in SMASH is 
shown by a hexagon. Finally, the five PSl discoveries (Lae 1, Tri H, Sgr H, Dra H, and Lae 3) are shown as stars. Bottom: The same for 
the size-Heliocentric distance plane. 
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Figure 5. Top: Particles of the fLaw fc Maiewskil l|2010bl'l N-body 
mod el of the Sgr stream, pro jected on the Sgr dSph co-ordinate sys¬ 
tem JMaiewski et al1l2003l) . Particles of the leading/trailing arm 
of the Sgr stream are shown in blue/green, whereas the body of 
the Sgr dSph is shown in orange. The position of Sgr II is repre¬ 
sented by the black star. Bottom: The same for the Heliocentric 
distance vs. Sgr dSph longitude plane. Sgr II clearly overlaps with 
the trailing arm. 



Heliocentric Distance (kpc) 

Figure 6. Heliocentric distance histogram of all particles in the 
ILaw &: MaiewsO ()2Q10bh model within 1 degree of Sgr IPs longi¬ 
tudinal position (no constraint was applied on the latitude). As in 
Figure [5] blue and green dots represent particles from the leading 
and trailing arm, respectively. The distance to Sgr H and its uncer¬ 
tainty are represented by the red star and the error bar and show 
that Sgr H is perfectly compatible with belonging to the trailing 
arm of the Sgr stream. 


discovered Laevens 1/Crater sys tem (r/,. = 20 ± 2 pc), 
recently confirmed to be a GC (|Kirby et al.ll201^ . It 
should however be noted that, recently, M31 satellite 
assu med to be GC hav e been discovered with similar 
sizes l|Huxor et al.ll201^ . although the nature of soin e 
of these systems is also ambiguous (IMackey et al.lIMT^ . 
The CMD of Figure 1 shows that the satellite hosts a 
clear blue HB, which is not a common feature of outer 
h a.lo GCs that tend to favor red HBs (see, e.g.. Figure 1 
of iLaevens et al.ll2014f l. Ultimately, spectroscopic follow¬ 
up and a measure of the satellite’s yelocity dispersion is 
necessary to fully confirm the nature of this satellite and 
whether it is dark-matter dominated. 

The location of Sgr II, ~ 15° from Sgr dSph and in the 
expected plane of the Sgr stellar stream is particularly 
interesting as it could point to an association. The fact 
that it lies 40-45 kpc behind Sgr dSph rules out a di¬ 
rect connectio n between the two satellite s but a compar¬ 
ison with the ILaw fc Maiewskil (l2010bll N-body model 
for the Sgr stream (Figure [5]) reyeals that Sgr II is lo¬ 
cated at the expected distance of model particles from 
the trailing arm of the Sgr stream stripped out of their 
host more than 3 Gyr ago. It is therefore likely that Sgr II 
was brought into the MW halo along with this part of 
the Sgr stream that has so far eluded detection, in a 
similar fashion to nume rous other MW halos GCs (e.g. 
ILaw fc Maiewskill20I0ar) . The fact that the sky location 
of Sgr II is slightly offset from this section of the model’s 
trailing arm is not necessarily surprising since a former 
Sgr dSph satellite is not expected to be as concentrated 
on the sky as its former stars in the model. In addition, 
the location of these older wraps of the Sgr stream is very 
poorly constrained in the model. In fact, the discovery of 
Sgr II and its association with the Sgr stream could add 
valuable constraints on the modeling of the Sgr stream 
once confirmed through radial velocities. 

longitu 

4.2. Draco II 

Draco II also has an ambiguous nature, although it is 
here driven mainly by the large uncertainties on its struc¬ 
tural parameters and distance, stemming from the faint 
nature of the object in the PSl data. With the current 
photometry, the satellite appears to share the properties 
of Kim 2 or Eridanus III, believed to be GCs (iKim et al.l 
l20I5bl : lBechtol et al.ll2015l : lKoposov et al.ll20I5ll . On the 
other hand, its uncertain properties are also completely 
compatible with thos e of Wil I, favored to be a DG 
(|Willman et al.l[2011ll . Here as well, deep photometry 
and/or spectroscopy are necessary to classify this sys¬ 
tem. 

Given the common connection between faint stel¬ 
lar systems and stellar streams, we investigate possi¬ 
ble associations of the satellite to known MW halo 
streams. The closest stream to Dra II is the GD-I stream 
(|Grillmair fc DionatosI [200^ . Placing the new satellite 
ont o the stream c oordin ate system (^i, defined 

by iKoposov et al.l (|2010ll . we find that it is located at 
(j)i ^ 17.1° and (j )2 ^ 3.8°. Though iKoposov et al.l (|2010[) 
do not have any measurements in this region (their mea- 

This coordinate system is a rotated spherical one, aligned 
with the stream’s coordinates, (pi and (p 2 represent longitude and 
latitude respectively. 
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surements range from (j>i = —60.00°, 4>2 = —0.64°±0.15° 
to (j)i = —2.00°, (j )2 = —0.87° ± 0.07°), the extrapolation 
of the orbit at the location of Dra II yield 4>2 ~ —2.7°, 
only 5-6° away from the satellite. However, the extrapo¬ 
lated distance of the stream reaches only ~ 12 kpc there, 
to be compared with Dra II’s ~ 20 kpc. Therefore, if 
the GD-1 stream does not significantly deviate from the 
iKoposov et al.l (I20I0I1 orbit, the current distance esti¬ 
mate for Dra II appears too high for a direct association. 

4.3. Laevens 3 

The small half-light radius of Lae 3 (7 ± 2 pc) places it 
well within the regime of GCs. With a relatively young 
age (~ 9 Gyr) and stellar populations that are not very 
metal-poor ([Fe/H] ~ —1.9), it would be natural to clas¬ 
sify Lae 3 as a “young o uter halo” GC found in the oute r 
region of the MW halo (|Mackev fc van den Berghll2nn5h . 
However, some caveats should be noted: the isochrone 
fit relies on the photometry currently available, which 
suffers from crowding. The presence of an RR Lyra star 
could be at odds with a young halo scenario since its 
presence would point to a system that is at least 10 Gyrs 
old. We find no possible connection of this new system 
with known stellar streams in the MW halo. 

5. CONCLUSION 

In this paper, we have presented the discovery of three 
new faint Milky Way satellites, discovered in the photo¬ 
metric catalog of the PSl 37r survey. The characteriza¬ 
tion of Lae 3 suggests that it is a GG, with properties 
similar to ‘young outer halo’ GGs. The two other sys¬ 
tems, Dra H and Sgr H, have an ambiguous classification. 
Dra H contains mainly main sequence stars, as well as a 
handful of probable RGB stars. It is very faint but its 
structural parameters are uncertain enough to prevent 
a classification as an extended GG or a compact dwarf 
galaxy. It is located close to the orbital path of the GD-1 
stream but its distance is in disagreement with the expec¬ 
tations of the stream’s orbit (~ 20 vs. ^ 12 kpc) and ap¬ 
pear to rule out an association. Finally, Sgr H is located 
in a part of the size-luminosity plane that contains no 
other known system, either more extended than known 
MW GCs, or more compact than known MW DGs in 
its luminosity range. Independently of its nature, Sgr H 
is particularly interesting as it lies at the expected lo¬ 
cation of the Sgr dSph stellar stream behind the bulge. 
In particular, the distance to the new satellite favors a 
connection with the currently undiscovered part of the 
trailing arm of the Sgr stream produced by stars stripped 
from the dwarf galaxy more than 3 Gyr ago. Ultimately, 
spectroscopic follow-up will be necessary to conclusively 
establish the nature of the last two satellites or confirm 
their connection with the GD-I and Sgr stellar streams. 
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